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K^,^, In an attempt to understand why the dominating magnetic field observed in the von-Karman-Sodium (VKS) dynamo experiment is 
axisymmetric, we investigate in the present paper the ability of mean field models to generate axisymmetric eigenmodes in cylindrical 
geometries. An a-effect is added to the induction equation and we identify reasonable and necessary properties of the a distribution so 
that axisymmetric eigenmodes are generated. The parametric study is done with two different simulation codes. We find that simple 
distributions of a-effect, either concentrated in the disk neighbourhood or occupying the bulk of the flow, require unrcalistically large 
values of the parameter a to explain the VKS observations. 
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^ 1 Introduction 

in ; 

Dynamo action generated by a flow of conducting fluid is the source of magnetic flelds in astrophysical 
' objects. Hoinogeno us dynamos also ha ve been observed in three laboratory experiments (Riga dynamo, 



Gailitis et a^.ll200d. Karlsruhe dyna mo, IStieglitz and Miilleiil200ll . Cadarache von-Karman-Sodium (VKS) 



experiment, Monchaux et al. 200?! ). The analysis of the dynamo effect beneflts from complementary ap- 



^ proaches of scientiflc computing and experimental studies. Indeed, experimental fluid dynamos offer an 
opportunity to test numerical tools which can then be applied to natural dynamos. While the flrst two 
experimental dynamos produced results in agreement with the predictions of numerical simulations the 
successful Cadarache von-Karman-Sodium experiment brought interesting unexpected features: (i) dynamo 
action is observed only with soft iron impellers and not with steel ones, (ii) the axisymmetric component 
dominates (i.e. the azimuthal mode (m = 0)) when the dynamo action occurs. 

Numerically simulating high permeability conductors embedded in conducting fluids is a challenging 
task and requires the development of new codes. We focus our attention in the present paper on a possible 
numerical answer to the second question: how can an axisymmetric magnetic field be generated? 

The mode (m = 1) was observed as predicted from kinematic dynamo s imulations in the Riga experiment 
using axisymmetric flows ( Gailitis and Freiberg 1980, Gailitis et al. 20041 ) and in the Karlsruhe experiment 
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using an anisotropic a-effect (iRadler et a/Jl200^) or simulations based on the realistic configuration con- 
siderering the 52 spin generators f.Tilgner,, 2002). Since num erical simulations prior to the VKS experiment 
(iMarie et fliJ 1200,1 lEavelet et aiJ 1200.4 IStefani et ffl/Jl200fil l also applied axisymmetric time averaged fluid 
flows, the mode (m = 1) was again predicted and utilized during the optimization process of the VKS 
impellers. (Recall that Cowling's theorem forbids the excitation of an axisymmetric magnetic mode from 
an axisymmetric velocity field.) 

Although the occurrence of the mode (m = 0) was a surprise, it does not contradict physics but rather 
demonstrates that the axisymmetric flow assumption is too simplistic. The counter-rotating impellers 
driving the VKS flow are responsible for a relatively high turbulence level compared to the flrst two dy- 
namo experiments, and a large spectrum of azimuthal mode s must be taken int o account. This may be 
done in various ways. One can use the nonlinear approach ( Bavliss et al\ 2007I ). however the Reynolds 
number which is achievable using direct numerical simulations is thousand times smaller than the ef- 
fective one, implying that turbulence is still poorly described. Alternatively, one can add some ad hoc 
non-axisymmetric flow modes in a kinematic code. Indeed, non-axisymmetric disturbances i n terni s of iii - 
ter mittent azimuthally drifting vorte x structures have been observed in water experiments by Marie ( 20031 ) 
and ide la Torre and Burguete ( 200?! ) but their influence on the dynamo process is unknown. A third pos- 
sibility is provided by the application of a mean fleld model where the unresolv ed non-axisymmetric small 
scale fluctuations are parametrized by an a-effect ( Krause and Radler 1980l ). A lthough the mean fleld 
approach is som ewhat controversial for large magnetic Reynolds numbers (see e.g. ICourvoisier et al.ll2006 



Sur et al. 20081 ). it is numerically far less demanding than direct numerical simulations and allows the ex- 



ploration of the parameter space more easily. This is the approach that we follow in the present study since 
it receives support from the second order correlation approximation (SOCA) as discussed in section ^ 13. 2i 
Assuming scale separation, the mean fleld approach parametrizes the induction action of unresolved 
small scale fluctuations via the a-effect. Our purpose is to determine whether an a-model can produce 
axisymmetric modes with realistic values of a. Using (•) to denote the averaging operation and primes for 
unresolved quantities, the most simple a-model states that 



{u' X b') = aB, 



where u' denotes the fluctuating velocity fleld, b' = /iqW the fluctuating magnetic flux density or induction, 
h' the fluctuating magnetic fleld, iiq the vacuum permeability, B the mean magnetic induction, and a the 
pseudo-tensor representing the a-effect. The dynamo action is analyzed by solving the kinematic mean 
fleld induction equation 



dB 
'dt 



V X (u X B + aB - r]V X B) 



(2) 



where u denotes a prescribed large scale velocity fleld and ij the magnetic diffusivity {ij = \/ ^qcj with a 
the electric conductivity). In the kinematic approach the velocity field is given and any back-reaction of 
the magnetic field on the flow by the Lorentz force is ignored. 

In case of a homogenous, isotropic, non-mirrorsymmetric small scale flow, the a-tensor becomes isotropic 
and reduces to a scalar. The a-effect provides an additional induction source by generating a mean current 
parallel to the large scale fleld. In general, the huge spatial extensions of astrophysical objects ensure that 
even a small and localized a-effect, p otentially supported by shear, is able to generate a large scale field 
(see for example ICharbonneaul l2005l for a review of current models of the solar magnetic cycle) . In the 
VKS experiment a possible source of the a-effect is the kinetic helicity caused by the shear between the 
outward driven fluid flow trapped betwee n adjacent impeller blades and the slower moving fluid in the 
bulk of the container (jPetrelis et a/.l 120071 ). Unfortunately, no helicity measurements are available for this 
experiment and the magnitude and the spatial distribution of the corresponding a-effect are not known. 
Although sign changes can occur wi thin the domain, th ere are two limits for the a distribution. Either it is 
concentrated in the impeller region (jLaguerre et a/.ll2008b .a). or it is uniformly distributed in the cylinder. 
We analyze these two extreme configurations in the paper. 

The present work attempts to identify essential properties of the a-effect that are necessary to generate an 
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axisymmetric magnetic field in VKS-like settings. In the VKS experiment, the critical magnetic Reynolds 
number is Rm^"^^* ~ 32 using soft iron impellers, where Rm is the magnetic Reynolds number defined as 
usual as the ratio between the stretching and the diffusive terms (see Eq. [6] below). Kinematic simulations 
are applied close to the onset of dynamo action and are used for the estimation of the magnitude of the 
a-effect necessary to produce the quoted dynamo threshold. 

One important numerical difficulty consists of implementing realistic boundary conditions on the mag- 
netic field between the conducting region and the vacuum. This is tackled in different ways in the two 
codes described in the Appendix. On e code is based on th e coupling of a boundary element method with 
a finite volume technique (FV/BEM, lOiesecke et al. 2008 ). the other is based on the co upling of a finite 



eleme nt method with a Fourier approximation in the azimuthal direction (SFEMaNS, iGuermond et al. 



20091 ). The two codes are validated by comparing the results they give on identical problems. 



The paper is organized as follows. In section [2l we consider a cylinder with no mean flow, the dynamo 
action is caused solely by an a^-effect. The simplified a-tensor is similar to the one used in the modeling of 
the Karlsruhe experiment. We compare periodic and non-periodic axial boundary conditions. We validate 
our two independent codes by observing that they give the same results up to non-essential approximation 
errors. In section[3l we model the geometry and the flow pattern of the VKS experiment. The mean velocity 
field is obtained from water experiments. The a-effect is investigated and the growth rates of the modes 
[m = 0) and [m = 1) are compared. Conclusions are proposed in sectional 



2 Axisymmetric a^-dynamos in cylinders 
2.1 OL-model 

We use cylindrical coordinates (r, 9, z) throughout the paper. We first focus our attention on a cylindrical 
mean field dynamo model with vanishing mean flow u = The a-effect is the only source of dynamo 
action and provides the necessary coupling between the poloidal and toroidal flelds which is essential to 
drive the dynamo instability. The induction equation ([2]) reduces to 

dB 

— = V X {aB -r]V X B). (3) 

To bound the parameter space, we consider a spatially constant but anisotropic a-effect that we assume 
to be nonzero only in the annulus Ri < r < Ro- We choose a simplified version of the a-effect modeling 
the forced helical fiow realized in the Karlsruhe experiment where a has the following tensor form: 

'ao 0\ 

a = I ao (4) 
0/ 

The motivation for this approach originates from the modeling of turbulent fiow structures in fast rotat- 
ing spheres and planetary bodies like the Earth. Simplified spherical a^-dynamos with scalar, isotropic, 
a-tensor and zero mean flow are known to exhibit properties that resemble essential characteristics of 
the Earth's magnetic field like dipole dominance and reversal sequences. However, more realistic spatial 
distributions of the a -effect like the anisotropic structure (jl]) usually generate equatorial dip ole dynamos 
( Giesecke et al. 2005 ) unless anisotropic turbulent diffusivity is accounted for ( Tilgner 20041 ) . 



In the rest of section [2] we analyze the dynamo action in a cylinder of radius R and vertical extension H. 
We evaluate the influence of the geometric constraints of the container (in terms of the aspect ratio H/R) 
and of the spatial distribution of the a-effect (homogenous or concentrated in the annulus Ri < r < Ro) 
on the dynamo threshold. We define the critical dynamo number 

^crit p 
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aspect ratio H/R 

Figure 1. Critical dynamo number in dependence of the aspect ratio H/R obtained from the FV/BEM algorithm. The transition between 
the modes (m = 0) and (m = 1) occurs at H/R si 0.79. 



where is the value of oq at which dynamo action occurs. 



2.2 Axisymmetric dynamos with a uniform a distribution 

The aspect ratio of the Karlsruhe device is H/R ~ 1 according to the overall size of the cylindrical container. 
When looking for periodic solutions for magnetic fields, since H would represent the half wavelength of a 
transverse dipole, we have to consider periodic cylinders of aspect ratio equal to two. In this case, the a- 
effect is uniformly distributed in the cylinder. Using the SFEMaNS code, we have ver ified that C^^^ — 



for th e mode (m = 1). This result agrees with the Karlsruhe analytical modeling from lAvalos-Zuhiga et al. 



The magnetic field observed in the Karlsruhe experiment was indeed perpendicular to the cylinder 



axis and described as an equatorial dipole. Note, that the critical dynamo number for the axisymmetric 
mode {C^^^ = 7.15) is larger than the one for the mode (m = 1) explaining that the mode (m = 0) was 
not observed. 

We now investigate whether the axisymmetric mode can become dominant when the aspect ratio H/R 
varies. We abandon the assumption of periodicity along the z-axis, and we consider the finite cylindric 
geometry. A uniform spatial distribution of a is still assumed. We compute the critical dynamo number 
for various aspect ratios in the range [0.2; 2.0]. The results obtained with the hybrid FV/BEM code are 
reported in Figure [H C omputations done w ith the SFEMaNS code give al most identical results. These 
results have been shown ( Stefani et al. 20091) to be simila r to those obtained by Avalos-Zuhiga et al\ ( 200?! ) 



applying the integral equation approach ( Xu et al. 20081 ) 



The main conclusion that we can draw from these computations is that the structure of the magnetic 
eigenmode passes from an equatorial dipole (m = 1) to an axial dipole (m = 0) as the container geometry 
passes from an elongated cylinder to a flat disk. The change of structure of the dominating mode is observed 
when the aspect ratio is about H/R ~ 0.79. If the aspect ratio of the Karlsruhe dynamo {H/R ~ 0.84 
when the curved ended pipes are disregarded) had been chosen 10% smaller, an axisymmetric mode might 
have been produced. This would have enhanced the analogy of the model with geomagnetism. Since the 
experiment is now dismantled, further numerical work on this point is rather academic. 

The eigenmodes associated with the modes (m = 0) for H/R = 0.5 and (m = 1) for H/R = 2 are shown 
in Fig. [21 The three panels on the left-hand side show the components of the axisymmetric magnetic 
field for the aspect ratio H/R = 0.5 and the six panels on the right-hand side show the eigenmode for 
H/R = 2. In the second case the geometric structure of the non-axisymmetric magnetic field is shown in 
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Radius Radius Radius Radius 

Figure 2. Geometric structure of the magnetic field of the a^-dynamo in a finite cylinder with uniform a-efTect (obtained from FV/BEM 
approach). Left side: H/R = 0.5, axisymmetric field at the marginal value C"'* = 9.8 (from top to bottom: Br, Bg, B^). Right side: 
H/R = 2, non-axisymmetric magnetic field in two meridional planes that differ by an angle of 90 degrees at C"'* = 5.3. (from left to 
right: Br, Bg, Bz). Solid curves denote positive values, dashed curves denote negative values. 



Table 1. Periodic and finite cylinders of aspect ratio H/R — 2 and different 
annular distributions of a. The value 7.15* results from unpublished data from 



Avalos-Zufiiea and PlunianI (120031. 


Case 


Periodic 


Periodic 


Finite 


Finite 


R^/R 


CS"t(^ ^ 0) 


CS"'(m = 1) 


CS"t(„ ^ 0) 


CS'-"(m = 1) 





7.15* 


4.8 


7.5 


5.2 


0.5 


9 


8 


9.4 


8.8 


0.8 


20 


20.2 


20.2 


20.4 



two orthogonal meridional planes. Note, that as the azimuthal wavenumber passes from (m = 1) (right) 
to (m = 0) (left), the typical scale in the axial direction passes from ~ 2 for H/R = 2 to less than 1 
for H/R = 0.5. When decreasing the aspect ratio, the geometric constraint on the axisymmetric current 
increases C^^^ by a factor of 2. 



2.3 Axisymmetric dynamos with an annular a. distribution 

We keep the aspect ratio H/R = 2 which produces a transverse field when a is uniform, and we show 
now that axisymmetric dynamos can be obtained when the distribution of a is annular, i.e. a{r) ^ 
if r E \Ri.-,R] and a(r) = if r G [0, i?j). This study is inspired by unpublished analytical results of 
Avalos-Zufiiga and PlunianI \2m± in periodic cylinders, where a change of structure of the dominating 



mode is observed for an annular distribution when i2j/i2 ~ 0.7. 

We compare the critical dynamos numbers obtained in periodic and finite cylinders for the same aspect 
ratio H/R = 2. We compare three cases: Ri/R = 0, 0.5, 0.8. The results are summarized in Table [H All 
the eigenmodes are steady at the threshold, the bifurcation is therefore of pitchfork type. Table [T] shows 
that, depending on the value of Ri/R, there is a change of structure of the critical mode. In the periodic 
case the mode (m = 0) is critical when Ri/R = 0.8 with C™* = 20 and the axial wavelength is A = 1. 
The mode [m = 1) is critical when Ri/R = 0.5 with C^"* = 8 and the axial wavelength is A = 2. The 
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(a) 





(a') 




Figure 3. Geometric structure of the magnetic field of the a^-dynamo in a periodic cylinder of aspect ratio H/R = 2 with an annular 
Q-effect: (a-c) (m = 0), C^"* = 20 and R^/R = 0.8; (a'-c') (m = 1), C^"' = 8 and R,/R = 0.5. Represented are the radial (a,a'), 
azimuthal (b,b') and axial (c,c') components in the plane 6 = 0. Data results from the SFEMaNS approach. 





(a) 



(b) 



(c) 




Figure 4. Geometric structure of the magnetic field of the o'^-dynamo in a finite cylinder of aspect ratio H/R = 2 with an annular 
a-effect: (a-c) (m = 0), C="* = 20.2 and Ri/R = 0.8 and for (a'-c') (m = 1), C^"* = 8.8 and Ri/R = 0.5 Represented are the radial (a, 
a'), azimuthal (b, b') and axial (c, c') components. Data results from the SFEMaNS approach. 



corresponding magnetic eigen-modes are shown in Figure El Note that, as expected, the magnetic field is 
concentrated in the region where a is nonzero. It is striking that the finite and the periodic results are 
similar when the ratio Ri/R is nonzero (see Table [U and Figure S]) . The only notable effect of the top 
and bottom boundary conditions is to concentrate the fields inside the box. For example, the structure 
of the mode {m = 0) shown in Figure [Ha-c) is concentrated around the equator and is dominated by a 
quadrupolar stationary state like a 's2t2' state. 

The main conclusion of this section is that axisymmetric a^-dynamos can be obtained in cylindrical 
geometries. In a finite cylinder with uniform distribution of a a dominant axisymmetric mode occurs when 
the aspect ratio is below the critical value H/R ~ 0.79. For H/R = 2, an annular distributi on of a can also 
generate an axisymmetric magnetic field if the a-effect is sufficiently localized. Note that iTilgnerl 
obtained also an axial dipole (in a sphere) using an anisotropic magnetic diffusivity {r]\\ > rj±). Compared 
to similar studies in spherical geometry, we have explored only a small region of the parameter space such 
as the spatial distribution of a, its symmetry properties, and its tensor structure. This study has validated 
our two different codes since they obtained identical results on simple configurations. This gives us some 
confidence before turning our attention to the VKS experiment. 
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Figure 5. VKS set-up design and mean-flow structure: (upper part) dimensions and technical details of the experimental set-up, with the 
copper vessel including the cooling system, the thin envelope separating the flow and the stagnant liquid sodium, the impeller with fltted 
blades and the shaft. Courtesy of the VKS team, (lower part) Simplified geometry used in the numerical simulations in non-dimensional 
units. The conductivity jump between the liquid sodium and the copper vessel is taken into account. The dashed line in the lower half 
denotes the region where dynamo action is supported by the localized «-effect (see section ^ 13.31 1. Note, that the simulations do not 
consider the region behind the impellers (lid layers). 



3 Mean field dynamos with a VKS type flow 

3.1 Experimental setup and mean (axisymmetric) velocity field 

In the von-Karman-Sodium experiment a flow of hquid sodium is driven by two counter-rotating impehers 
located at the opposite ends of a cylindrical vessel. Self-generation of a magnetic field occurs if the magnetic 
Reynolds number exceeds the critical value Rm'^'^^* ~ 32. It is important to note though that the dynamo 
is activated only when impellers made of soft iron are used. A sketch of the experimental setup is shown 
in Figure [5l The vessel is filled with liquid sodium. Two counter-rotating impellers are located at the 
endplates of the vessel each fitted with eight bended blades. The flow generated by the rotating blades is 
of von-Karman type and has been extensively investigated in Ravelet et al. ( 20051 ). This so called s2t2-flow 
basically consists of two poloidal and two toroidal cells. 

Throughout this paper, units are non-dimensionalized by setting the radius of the flow domain to R = 1 
(corresponding to 0.206 m in the experiment, see Fig. [5|) and the magnetic diffusivity to rj = 1 (correspond- 
ing to 7/Na ~ 0.1 m^/s in the experiment). Dimensional values for the (maximum of the) velocity or the 
magnitude of the a-effect are then obtained in units of ms~^ by a multiplication with a factor close to 0.5. 
To facilitate the dynamo action the flow region is surrounded by a layer of stagnant sodium of thickness 
0.4 which is enclosed by a solid wall of copper of thickness 0.2. The fluid velocity is zero when r > R = 1. 
The copper wall is modeled by a high conductivity region, i.e. ficu ~ 4.5(TNa when 1.4 < r < 1.6. The flow 
is characterized by the magnetic Reynolds number 

Rm = ^^21^ (6) 
V 

where ?7max is the maximum of the fluid velocity. A typical example for a flow field realized in a water 
experiment in the VKS geometry is shown in Figure [H It is time-averaged, axisymmetric and symmetrized 
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Figure 6. Axisymmetric velocity field obtained from water experiments performed by iRavelet et al\ l|2005l ) and which is applied in 
present kinematic simulations. Grey scale coded structures represent the azimuthal velocity field ug and the arrows represent the 
poloidal components {ur,Uz). 




Figure 7. 
Bf,. 



Iso-surface of the magnetic energy density at 40% of the maximum value. The grey scale represents the azimuthal component 



with respect to the equator. We utiUze this mean velocity field in the kinematic induction equation ([2]). 

Using our two codes, we find that the critical Reynolds number for the mode (m = 1) for this setting is 
approximately Rm'^"* ~ 45 without any a-effect. The geometric structure of the eigenmode obtained for 
Rm = 60 is shown in Figure[71 The picture shows the iso-surface of the magnetic energy density (2/io)~^ 
at 40% of the maximum value. The grey scaled map on the iso-surface represents the azimuthal component 
of the magnetic induction. This map shows that the mode (m = 1) dominates. The well known embracing 
banana- like structure is the dominating feature in the central region. The accumulation of magnetic energy 
close to the equator at the interface between the flow domain and the surrounding stagnant layer of fluid 
is another striking feature. Because of the equatorial symmetrization o f the velocity field , the resulting 
eigenmode does not exhibit any azimuthal drift and remains stationary ( Marie et al. 20031 ). 
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Figure 8. Sketch of the «-model in the VKS experiment: between two impeller blades, (left) the fluid velocity (thin arrow and outwards 
symbol) acting on a given azimuthal mean magnetic field B (thick arrow) results in a secondary magnetic field 5 (middle) and generates 
an electromotive force u X b (right). 



3.2 a-model 

We now want to explore the implications of the suggestions that the strong shear between the flow trapped 
between the impeller blades and the slower rotating fluid neighboring the impeller results in a helical 
outward driven flow. The induction efl'ects of such a flow can be parametrized by an a-effect as illustrated 
in Figure [HI In this sketch an azimuthal electromotive force is produced from the cross-product of the 
flow u between the blades and the magnetic field b obtained from the distortion of an applied azimuthal 
magnetic field. The resulting electromotive force is parallel to the applied magnetic field but of opposite 
direction. This corresponds to a negative coefficient agg in the a-tensor. Applying a vertical magnetic field 
leads to a vertical electromotive force again parallel to the applied magnetic field and of opposite direction. 
This corresponds to a negative coefficient a^z in the a-tensor. 

In order to estimate the value of the coefficients agg, azz-, we first estimate the magnetic Reynolds 
number of the fluid flow trapped between the impeller blades. Let C/disk be the rotation speed at the rim 
of the disc. The maximum fluid velocity is obtained just below the upper im peller and has been , measured 



to be C/max ~ i^f^disk with V = 0.6 denoting an impeller efficiency parameter ([Ravelet et a/.ll2005l ). Then an 
estimate of the recirculation flow intensity is u ~ (f^disk — t^max)/2. This leads to u/C/disk ~ (1 — 1^)/2 = 0.2. 
The height of a blade is /i i^/20 « (1.8/20)i? = 0.09ii, R being the disc radius and H the distance 
between the discs. Reminding the deflnition of the global magnetic Reynolds number from Eq. ([6]), Rm = 
UmaxR/v ~ t^UdiskR/i] we define a local magnetic Reynolds number that characterizes the fiow between the 
blades as Rmbiadc ~ uh/rj. The combination of both definitions yields Rmbiade ~ iu/Udisk){h/R)(Rui/i') 
leading to Rmbiade ~ 1 for Rm ~ 32. Therefore it seems reasonable to apply a low magnetic Reynolds 
number approxi mation to estimate the v alue of the a-effect in the experiment according to the SOCA 
approximation ( Krause and Radler 1980l ). This gives aeg ~ a^z ~ uurh/rj. Assuming that the radial 



outward flow Ur is of the same order as u, for ry = 0.1 m^s ^ ar id = 0.2 m w e obtain age ~ Uzz ^ 5 ms ^ . 
This estimate exceeds the value (a ~ 1.8 ms~ ) reported by lLaguerre et al. (l2008bl ) where the efficiency 



factor V was not considered in the definition of Rm. Note, that this is only an evaluation of the order of 
magnitude and, as the fiow is highly turbulent, the optimum value of a estimated above probably is not 
realized in the experiment. 

We henceforth focus our attention on agg. It is indeed this coefficient that generates the poloidal com- 
ponent of an axisymmetric magnetic field from its toroidal component. In turn, the toroidal component 
is efficiently generated by twisting and stretching the poloidal field component via a shear flow (w-effect). 
The impact of an Uzz coefficient will be discussed at the end of the next section. 
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Figure 10. Field amplitude growthrates for the localized a-effect. 
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3.3 Localized otge-effect results 

As the spatial distribution of a is not known a priori, we start by performing simulations with the a-effect 
restricted to the impeller region. More specifically we set 



a ( 1 + 




tanh 



2-Ztor 



tanh 



if < r < i? 
if i? < r 



(7) 



where a can be positive or negative, Az = 0.05, ztop = Zbot = iO.8, and i? = 1 is the radial extension of 
the impeller region. Equation ([7]) specifies a smooth transition of ago between a maximum value a in the 
impeller region and a vanishing value in the bulk of the container. We henceforth refer to this model as 
the localized a-effect. The axial profile of agg{z) is shown in Figure [H 

Figure [10] shows the growthrates of the magnetic field as a function of a for various magnetic Reynolds 
numbers. The solid curves represent the growthrate for the axisymmetric mode {m = 0) and the dashed 
curves represent the results for the mode (m = 1). We observe two distinct dynamo regimes that differ 
by the resulting field geometry. The main characteristics of the dynamo solutions can be summarized as 
follows. For a < —50 the eigenmode is axisymmetric whereas for a > — 15 the mode (m = 1) dominates. 
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Figure 12. Critical magnetic Reynolds number in dependence of the magnitude of the o-effect. The black curve represents the results 
obtained from the SFEMaNS scheme and the grey curve denotes the results of the FV/BEM approach. The left part (solid curves) shows 
pj^j^crit £qj. ^jjg axisymmetric mode (m = 0) whereas the left part (dashed curves) is obtained for a dominating mode (m = 1). 



We observe in Figure [TOl that the growth rates of both modes are negative in the range [—50, —10]. This 
interval might become smaller as Rm increases, but it seems that there exists a non-empty range of a for 
which no dynamo action is possible at all. For instance, at a = —20 no dynamo is obtained for Rm up to 
300 and the corresponding growthrates even decrease as Rm increases (see Figure [TTI) . 

The existence of these two regimes becomes more obvious on Figure [12] where we show the critical 
magnetic Reynolds number as a function of a. (Rm'^'^'* is estimated by linear interpolation of two adjacent 
growthrates obtained close to the dynamo threshold.) The black curves represent the results provided 
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by the SFEMaNS code and the grey curves represent the results from the FV/BEM code. Slight but 
systematic deviations between both codes are observed especially for large negative values of a. Most 
probably the disagreement is the result of a small difference in the localization of a and/or the velocity 
field that is caused by the staggered mesh definition in the FV/BEM scheme. The present work does not 
intend to perform a benchmark examination of different numerical schemes so that the deviations are not 
investigated in detail as the general trend and the critical values of the magnetic Reynolds number are in 
rather good agreement. 

By examination of the left branch of the graph in Figure [12] we estimate that agg ~ —115 is necessary to 
obtain an axisymmetric dynamo at Rm ~ 32. In physical units this value corresponds to a ~ 57.5 ms~^, 
which exceeds the estimated magnitude of a by more than a factor of 10. Furthermore, this value is also 
four times larger than the maximum fluid velocity in the bulk of the cylinder. It is hard to believe that 
such a large a-effect is realized in the experiment. 

We performed additional simulations with an Oizz coefficient located at the impellers. The o^-tensor is 
therefore given by (0,06)61,022) with a tanh localized function (see Eq. [Tj) such that aQQ = azz- We have 
found, for Rm = 60, Oc = -42.60 with SFEMaNS and -44.6 with FV/BEM. These values are to be 
compared with Uc = —83.40 for Rm = 60 obtained with a localized (0, a^g, 0) a-tensor. This reduction by 
a factor 2 is not enough for approaching realistic values of a. 

3.4 Uniform cx-effect results 

In this section we assume that the a-effect is uniformly distributed in the flow region (0 < r < 1, —0.9 < 
z < 0.9). Here, only the age component is assumed to be nonzero. 

Figure [13] shows the growthrates as a function of a for Rm = 30, 50 and 100. These results have been 
obtained with the FV/BEM code. We have verified on a few simulations (not reported here) that the 
SFEMaNS code produces similar results. 

We see that for large negative values of a the mode (m = 1) is suppressed in all the cases and the 
corresponding growthrates remain negative. The growthrate of the axisymmetric mode becomes positive for 
sufficiently large magnetic Reynolds numbers. Contrary to what is observed when the a-effect is localized, 
positive growthrates of the axisymmetric mode are obtained for positive values of a. For a ^ the 
growthrate becomes very sensitive to small changes in a, and a small interval exists around a ~ +5 where 
the critical magnetic Reynolds number of the modes (m = 0) and (m = 1) are rather close together. This 



May 25, 2009 7:19 Geophysical and Astrophysical Fluid Dynamics gicscckc'AxisymMFDynVKS 

13 

200 



150 



o 

^ 100 



50 





-100 -80 -60 -40 -20 20 40 

a 

Figure 14. Rm*^"^'' for different o-distributions. Black curves: localized a-effect. Grey curves: homogenous o-effect. Solid lines denote the 
results for the mode (m = 0), dashed lines denotes the results for m = 1. 

might be a promising possibility to obtain an axisymmetric field with more reasonable values for \a\. This 
alternative needs to be examined in more details because positive values for a can be obtained in the 
equatorial layer where large scale intermittent vortices have been observed ( 

For —15 ^ a ^ —7, again, there exists a region that is characterized by an extremely high criti- 
cal magnetic Reynolds number, if a dynamo could exist at all. Figure [TH shows the critical magnetic 
Reynolds as a function of a for both the localized and the uniform a-effect. The behavior of Rm'^"* for 
either a-distributions exhibits similarities like a strong tendency to diverge around a certain value of a. 
Approximating the curves with a power law Rm'^'^'* oc (a — ag)*^ gives a diverging behavior around Og ~ —18 
for the localized a-effect and a^ ~ —7.4 for the homogenous distribution. 

The critical range without dynamo action becomes significantly smaller when the a-effect is assumed to 
be uniform. From the negative branch (a ^ 0, left graphs in Figure [T4l) we see that when the a distribution 
is uniform, the value of |a| which is necessary to produce an axial dipole field at a certain Rm'^"* is 
significantly lower than when the a-effect is localized. To obtain an axisymmetric field at Rm'^'^'* ~ 32 we 
need to set a « —115 when the distribution is localized and a ~ —35 when the distribution is uniform. 
Comparing the volume fractions occupied by the a-effect in both configurations (which differ by a factor 
of about 10), it becomes clear that the a-effect within the impeller region operates more effectively than in 
the remaining active zone. However, even when the a-effect is assumed to uniformly penetrate the entire 
fiow domain the magnitude of a that is required to obtain the mode (m = 0) at Rm'^'''* ~ 32 is still 
|a| |C/max| ~ 16 ms^^ which remains unrealistic high. 




Mariell2003l l. 



4 Conclusion 

The present study is an illustration in the context of the fiuid dynamo problem of the interaction between 
numerical and e xperimental approaches. An account of previous stages concerning the VKS experiment 
may be found in Leorat and Nor j ( 20081 ). We have shown that in the framework of a mean-field model, 
axisymmetric a^-dynamos are possible in cylinders embedded in vacuum. It seems always possible to find 
a configuration (either determined from geometry or by a certain spatial distribution of a) which is able 
to generate an axisymmetric magnetic field. 

In a VKS-like configuration, the combination of an axisymmetric fiow and an a-effect can produce 
axisymmetric magnetic modes as well. Our simulations with simple profiles of a point out however that 
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large and unrealistic values of a are necessary to explain the VKS experimental results. We could think of 
more complex a distributions with negative values between the blades and positive values in the equatorial 
layer. This is left for future work since a realistic assessment would require to measure the kinetic helicity 
distribution in a water model of the VKS device. 

We note also that, even if realistic values of the magnitude of the a-effect had been obtained, one should 
have then to explain the non-existence of the dynamo action when using steel impellers. Indeed, at first 
sight, steel blades produce a flow, including kinetic helicity, close to the one driven by the soft iron blades. 
The role of the ferromagnetic material to obtain the dynamo action appears to be a critical issue and 
deserves further experimental and numerical investigations. 
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Appendix ANumerical methods 
A.l SFEMaNS algorithm 

The SFEMaNS acronym stands for Spectral Finite Element method for Maxwell and Navier-Stokes equa- 
tions. This code is designed to solve the MHD equations in axisymmetric domains in three space dimensions. 
To simplify the presentation we restrict ourselves in this Annex to the kinematic induction equation. 

Let us consider a bounded domain $7 C with boundary F = The domain i7 is assumed to 
be partitioned into a conducting region (subscript c) and an insulating region (subscript v): Q = 
Q^, ricHO^ = 0. The interface between the conducting region and the nonconducting region is denoted by 
S = dilc n dQv To easily refer to boundary conditions, we introduce Fc = F n dQc, F^, = F n dQy Note 
that F = F„ U Fc. 



A. 1.1 The PDE setting. The electromagnetic field in the entire domain i7 = U 0^ is modeled by the 
Maxwell equations in the MHD limit. 



' txdtR = -VxE, in VL 
VxH -- 



(t(E + ■UX//H) + j**, inVLc 
0, in fi, 



(A.l) 



VE = 0, in J7„ 

E X n|r = a, ^\t=Q 



Ho, in Q.C 



where n is the outward unit normal on F. The independent variables are space and time. The dependent 
variables are the magnetic field, H = B//i, and the electric field, E. The data are the initial condition. 
Ho, the boundary data, a, and the externally imposed distribution of current, j*. The data are assumed to 
satisfy all the usual compatibility conditions, i.e. V(^Ho) = 0. The physical parameters are the magnetic 
permeability, and the conductivity, a. For the sake of generality, the permeabilities in each domain 
{pi'^,^'") are distinguished, but we take fi'^ = fi" = fiQ in the applications considered in the paper. 
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A. 1.2 Introduction of 4> and elimination of E. We assume that the initial data Hq is such that 
VxHo|n„ = 0. We also assume that either is simply connected or there is some mechanism that ensures 
that the circulation of H along any path in the insulating media is zero. The condition VxH[q^ =0 then 
implies that there is a scalar potential (p, defined up to an arbitrary constant, such that H|q^_ = V0. 
Moreover, we can also define (po such that Ho|n„ = V</>o. We now define 



H 



in 

V(j) in Q^, 



/i*^ in 
/i^ in fit,. 



(A.2) 



and we denote n'^ and the outward normal on dQc and dQ^, respectively. It is possible to eliminate the 
electric field from the problem and we finally obtain: 



H^dtAc/) = 


in Qy 


{R^^a-\V X - j") - u X ^'^H^) x = a 


on Tc 


fi^dn^{dt(p) = -n^-Vx(n''xa), 


on Ty 


H^xn^ + Vc/.xn" = 


on S 


//'=H^-n^ + ^^V</.-n^ = 


on S 


Il%=o = Hg, 0|t=o = ^0- 





(A.3) 



A.1.3 Weak formu lation. The weak formulation of (IA.3P that we want to use has been derived in 
Guermond et all ( 200?! ). We introduce the following spaces: 



L 

X 



{(b,v9) GL2(fi,)x/7i„(fi,)}, 

{(b,v9) eHeuri(J^c)xF^i„(fi.); (b x + X n'^) |s = 0}, 



(A.4) 
(A.5) 



and we equip L and X with the norm of Ij'^{Q,c)^H^{Q,y) and Hcmii^c) ^ i^v) , respectively, //^^^(fi^) 
is the subspace of H^{Clv) composed of the functions of zero mean value. The space Hcuri(^^c) is com- 
posed of the vector-valued functions on fic that are componentwise L^-integrable and whose curl is also 
componentwise L^-integrable. 

We are now able to formulate the problem as follows: Seek the pair (H'^,</)) G L^((0, +00); X) n 
L°°((0, +00); L) (with dtil^ and dt(j) in appropriate spaces) such that for all (b, 95) G X and t E (0,+cx3), 

' H^|i=o = Hg; V0|t=o = 

/ [^^(atH^)-b + ((i?ma)-i(VxH^-j*') -itx^^H^)-Vxb] + / fi''{dtV(l))-Vip 

r , (A.6) 

+ J ((i?mf7)"^(VxH^-j^) -wx/i^H^)-(bxn^ + Vv^xn") 

= J (axn)-(bxn) + J (axn)-(V93Xn). 

The interface integral over T, is zero since bxn'^ + Vipxn" = 0, but we nevertheless retain it since it 
does not vanish when we construct the nonco nforming finite element approximation in ^A.1.4[ 

It has been shown in Guermond et al. ( 20071 ) that ()A.6P is equivalent to ()A.3p . Observe that the boundary 
conditions on Ty and in (jA.Sp are enforced naturally in ()A.6p . The interface continuity condition 
H'^xn'^ + V^xn'' = is an essential condition, i.e. it is enforced in the space X, see ()A.5p . One originality 
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of the approximation technique introduced in iGuermond et o.l\ (jioO^) and recahed in §A.1.4I is to make 
this condition natural by using an interior penalty technique. 



A. 1.4 Finite element approximation. We approximate ()A.6P by means of finite elements in the merid- 
ian section and Fourier expansions in the azimuthal direction. 
The generic form of approximations of H*^ and (j) is 

M 

f{r,e,z,t)= tt{r,z,t)e^', fjl{r,z,t)=f^''{r,z,t),ykeO;M, (A.7) 

k=-M 



where i^ = —1 and M + 1 is the maximum number of complex Fourier modes. The coefficients /^(r, z,t) 
take values in finite element spaces. We use quadratic Lagrange finite elements to approximate the Fourier 
components of H*^, i.e. the three components of the magnetic field, {H^, Hg, H^), are continuous across 
the finite elements cells and are piecewise quadratic. The approximation space for the magnetic field is 
denoted X^. Similarly, the Fourier components of the magnetic potential are approximated with quadratic 

Lagrange finite elements. The approximation space for the magnetic potential is denoted X^. No continuity 
constraint is enforced between members of X^ and X^. 



Th e Maxwell equation is approximated by using the technique introduced in IGuermond et al\ (|2007l . 
20091 ). The main feature is that the method is non-conforming, i.e. the continuity constraint (bxn'^ + 
\/(pxh")\Y; = in X (see (]A.5|) ) is relaxed and enforced by means of an interior penalty method. 

We use the second-order Backward Difference Formula (BDF2) to approximate the time derivatives. The 
nonlinear terms are made explicit and approximated using second-order extrapolation in time. Let At be 
the time step and set t"" := nAt, n > 0. The solution to the Maxwell equation is computed by solving for 
■jj^c,n+i -j^ -j^H g^^^ (^"+1 in xf^ SO that the following holds for all b in X^ and all ip in xf^ 



b + (i?mfT)"^VxH'='"+i-Vxb 



At 



{- [ ((i?„iCT)"^(VxH^'"+i -j^) -ux/^^H*)-(bxn^ + Vt^xn") 
+ <7((H^'"+\ 0"+^), (b, ip)) + .(H^'«+\ b) 
[ {uxfi^U* + {R^a)-^j')-Vxh+ [ (axn)-(bxn)+ [ (axn)-(V(/?xn), (A.^ 



where we have set DH^'"+i := i(3H^'"+^ - 4H'^'" + H^'""!), := i(3(/)"+i - 4(/)" + (^""i), and 

5((H^'"+\ (b,^)) := [3 J2 I (H^'"+^xn^ + V0"+^xn^)-(bxn'^ + V99xn''), (A.9) 

s(H^'"+\b) :=7 /" V(^"H'='"+^)V(/x^b). (A.IO) 



The purpose of the bilinear form g is to penalize the quantity 11"^'"+ xn'^ + V^"'^ xn'' across S so that it 
goes to zero when the mesh-size goes to zero. The coefficient (3 is user-dependent. We usually take (3=1. 
The purpose of the bilinear form s is to have a control on the divergence of H*^. 
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A. 2 3D finite-volume/boundary-element-method (FV/BEM) 



A. 2.1 Finite volume method. A finite volume (FV) approach provides a robust grid based scheme for 
the solution of the kinematic induction equation. The local spatial discretization on a regular grid is easy 
to implement and allows a fast numerical solution of the induction equation in three dimensions. Physical 
quantities are defined at distinct locations on grid cells that are obtained from a regular subdivision of the 
computational domain. The components of the at a grid ceh labeled by (ix, iy, iz) 

are defined at the center of the cell faces and the values are interpreted as the average of the magnetic 
field on the specific cell face. The field update at a time step n + 1 requires the discretization of Faraday's 



law d+B 



-V X E. This implies the computation of the electric field E. 



which is defined on the 



edges of a grid cell so that the localizations of the components of E are slightly displaced with regard 
to the components of B. Additional computational efforts occur as the computation of E requires the 
reconstruction of u and B on the edge of a grid cell. For moderate magnetic Reynolds numbers it is 
sufficient to apply a simple arithmetic average to interpolate u and B on th e edge. For larger magnetic 
Reyn olds numbers, however, more elaborate schemes have to be applied (e.g. Ziegler 2004I . Teyssier et al. 

The complications that arise by the definition of a second, staggered mesh are essentially outweighed 
by the maintenance of the divergence free condition and the conservation of the fiuxes across interfaces 
between neighboring cells which are intrinsic properties of the specific finite volume approach. 

To relax the constraints of the time step restriction an implicit solver for the diffusive part — V x (ryV x B) 
of the induction equation is applied. The full scheme for the semi-impli cit field update at tim e step n + 1 
is second-order in time and is summarized by the following expression ( Keppens et al. 19991 ): 



B 



n+l 



B^ + AtF'^P 



2 ^ 



(A.ll) 



Here F'^^^ denotes the discretized operator accounting for the terms of the induction equation that have 
been made explicit {(x u x B and aB), F™p denotes the discretized operator accounting for terms made 
imphcit (oc ??V x B) and F = F™p + F'^^p. 



A. 2. 2 Insulating boundary conditions. Insulating domains are characterized by a vanishing current 
J oc V X S = so that, assuming that the vacuum domain is simply connected, B can be expressed as the 
gradient of a scalar magnetic potential B = — V<^. The potential is a solution to the Laplace equation 



A<? = 0, ^ ^ 0(r"2) for r oo. 



(A.12) 



The computation of B at the boundary requires the integration of = 0. An effective approach 

to compute the potential ^ and the corresponding boundary field for insulating boundary conditions is 
provid ed b y the boundary element m ethod (BEM). This procedure has been proposed in llskakov et al. 
tOOi) andllskakov and Dormvl (120051 ). and was recently modified and applied to various dynamo problems 
riesecke et al\ 20081 ). We now give a short description of the technique. 

Consider a volume that is bounded by the surface F and let d/dn = n • V denote the outward 



normal derivative. After applying Green's second theorem including some straightforward manipulations, 
the magnetic potential ^ is shown to satisfy the following integral equation: 



1 



(A.13) 



where G{r,r') is the Green's function, or fundamental solution, which fulfills AG{r,r') = —6{r — r') 
and is explicitly given by G{r,r') = — (47r \r — r'\)^^. (9„$ = —B^ yields the normal component of B on 
dT which is known from the finite volume method as described in the previous paragraph. The tangent 
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components of the magnetic field at the boundary are computed by: 

B^ = er-B= -er-V<P{r) = 2jer-(^{r')\/r^^^^^+B'\r')V/^{r,r')Jdr{r') (A.14) 

r 

where e-j- represents a tangent unit vector on the surface element dT(r'). It can been shown that Eqs. ()A.13P 
and (|A.14p contain all the contributions if all the field sources are located within the volume enclosed by 

r. 

The discretization of Eq. (jA.lSh and Eq. ()A.14p leads to an algebraic system of equations which finally 
determines B^ at a single point on the surface by performing a matrix multiplication that connects the 
normal components B^ at every surface grid cell of the computational domain: 

N 

BJ = Y,M^JB]. (A.15) 

i=l 

Eq. (jA.lSP introduces a global ordering of the quantities B'^ and i?" defined by an explicit mapping of 
the boundary grid-cell indices {ix, iy, iz) on a global index i = 0, 1, 2, • • • , where N = 2 ■ [nz ■ ny + nz ■ 
nx + ny ■ nx) represents the total number of boundary elements. The matrix elements Mij are computed 
numerically applying a standard 2D-Gauss-Legendre Quadrature method. In general the matrix M. is 
dense and requires a large amount of computational resources. However, M. only depends on the geometry 
of the problem and therefore has only to be computed once. 
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